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Abstract. We present a survey of the known cosmological and black hole solutions 
in ghost-free bigravity and massive gravity theories. These can be divided into three 
classes. First, there are solutions with proportional metrics, which are the same as in 

£L| General Relativity with a cosmological term, which can be positive, negative or zero. 

Secondly, for spherically symmetric systems, there are solutions with non-bidiagonal 
metrics. The g-metric fulfills Einstein equations with a positive cosmological term and 
a matter source, while the f-metric is anti-de Sitter. The third class contains solutions 

^. | with bidiagonal metrics, and these can be quite complex. The time-dependent solutions 

00 ' describe homogeneous (isotropic or anisotropic) cosmologies which show a late-time 

self- acceleration or other types of behavior. The static solutions describe black holes 
with a massive graviton hair, and also globally regular lumps of energy. None of these 
are asymptotically flat. Including a matter source gives rise to asymptotically flat 



f»^ \ solutions which exhibit the Vainshtein mechanism of recovery of General Relativity in 

(*C) ■ a finite region. 
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1. Introduction 

The recent discovery of a multiparameter class of ghost-free massive gravity theories 
by de Rham, Gabadadze and Tolley (dRGT) [TJ, and of its bigravity generalization by 
Hassan and Rosen (HR) [2], has revived the old idea that gravitons can have a small 
mass m [3] (see jl] for a recent review). Theories with massive gravitons were, for a 
long time, considered as pathological, mainly because they exhibit the Boulware-Deser 
(BD) ghost - an unphysical negative norm state in the spectrum [5]. However, there 
are serious reasons to believe [6], El [HI El QUI ED] that the theories of [TJ [2] are free of 
this pathology. Hence, they can be considered as healthy models, and it is interesting 
to explore their physical consequences. 

Theories with massive gravitons have been used in order to explain the observed 
acceleration of our universe [321 EG]. This acceleration could be accounted for by 
introducing a cosmological term in Einstein equations, however, this would pose the 
problem of explaining the origin and value of this term. An alternative possibility is to 
consider modifications of General Relativity (GR), and theories with massive gravitons 
are natural candidates for this, since the graviton mass can effectively manifest itself as a 
small cosmological term [H] . This motivates studying cosmological and other solutions 
with massive gravitons. 

Theories with massive gravitons are described by two metrics, g^ v and f^. In 
massive gravity theories the f-metric is non-dynamical and is usually chosen to be flat, 
although other choices are also possible, while the dynamical g-metric describes massive 
gravitons. In bigravity theories both metrics are dynamical and describe together two 
gravitons, one of which is massive and the other one is massless. The theory contains 
two gravitational couplings, n g and Kf, and in the limit where /t/ vanishes, the f- 
equations decouple from g, so that the f-metric is determined by its own dynamics. 
The g-equations contain f, and the g-theory then can be viewed as a massive gravity 
with a fixed background metric f. Therefore, massive gravity theory is a special case of 
bigravity, and so when one studies bigravity solutions one finds, in particular, all massive 
gravity ones. One should emphasize, however, that in the Kf — > limit the f-metric 
does not necessarily become flat and the bigravity solutions do not always reduce in this 
limit to solutions of the standard dRGT massive gravity with flat f. 

The known bigravity solutions can be divided into three types. First, there are 
solutions for which the two metrics are proportional in the same coordinate system: 
fiiv = C 2 g^ v . This is only possible if C fulfills an algebraic equation with constant 
coefficients, and if the matter sources in the two sectors are fine-tuned to be proportional 
to each other. The g-equations then reduce to the standard Einstein equations with a 
matter source and a cosmological term A g ~ m 2 , which can be positive, negative or zero. 
Therefore, one obtains in this way all GR solutions, in particular, for C = 1, those of 
vacuum GR. For A g > there are cosmological solutions which approach the de Sitter 
space at late times and describe the cosmic self-acceleration. On the other hand, none 
of the proportional solutions, apart from the trivial one with f^ u = g^ u = i]^ v , fulfill 
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equations of the massive gravity theory with flat f. 

Secondly, imposing spherical symmetry, there are solutions described by two metrics 
which are not bidiagonal, although they can be made separately diagonal when expressed 
in two different frames. There is a nontrivial consistency condition relating these two 
frames, and as soon as it is fulfilled, the two metrics formally decouple one from the 
other and each of them fulfills its own set of Einstein equations with its own cosmological 
and matter terms. In some cases, the f-metric can be chosen to be anti-de Sitter, with 
Af ~ kJ, hence it becomes flat when Kf — y 0, and so the dRGT massive gravity is 
naturally recovered. 

These solutions exist both in the dRGT massive gravity and in the HR bigravity 
theories. They describe all known massive gravity black holes, in which case the g-metric 
is Schwarzschild-de Sitter. They also describe all known massive gravity cosmologies, 
for which the g-metric is of any (open, closed and flat) standard Friedmann-Lemaitre- 
Robertson- Walker (FLRW) type. For these solutions the ordinary matter dominates 
at early times, when the universe is small, while later the effective cosmological term 
A g ~ m 2 becomes dominant, leading to a self-acceleration. Such solutions had been 
first obtained without matter, in which case they describe the pure de Sitter space 
[T5| IT6] . Later the matter term was included, [T7J HHJ [191 I2H], first only for special 
values of the theory parameters, and then in the general case [2TI 122] . If the massive 
gravity theory were indeed the correct theory of gravity, these solutions could describe 
our universe. However, perturbations around these FLRW backgrounds are expected to 
be inhomogeneous - due to the non-diagonal metric components, although this effect 
should be suppressed by the smallness of m [18]. For this reason solutions of this type are 
sometimes called in the literature 'inhomogeneous', or solutions with 'inhomogeneous 
Stiickelberg fields'. 

For solutions of the two types discussed above the metrics are the same as in 
GR and the graviton mass manifests itself only as an effective cosmological term. By 
contrast, the third type of the known bigravity solutions, where the two metrics are 
bidiagonal in the same frame but not proportional, leads to more complex equations 
of motion which usually require a numerical analysis. The FLRW solutions can be 
self-accelerating, but can also show more complex behaviors [23l |24[ |25| [26] . There are 
also anisotropic cosmological solutions, and these show that the generic state to which 
the universe approaches at late times is anisotropic, and the anisotropy energy behaves 
similarly to the energy of a non-relativistic matter [27J. The static vacuum solutions 
with bidiagonal metrics can be rather complex and describe black holes with a massive 
graviton hair and also gravitating lumps of energy [28J. None of these solutions are 
asymptotically flat. The asymptotic flatness is achieved by including a regular matter 
source, which gives rise to solutions which exhibit the Vainshtein mechanism of recovery 
of General Relativity in a finite region of space [29] . 

Below, a more detailed description of the currently known bigravity and massive 
gravity solutions is given. 



M.S.Volkov 4 

2. Ghost-free bigravity theory 

The theory of the ghost-free bigravity [2] is defined on a four-dimensional spacetime 
manifold equipped with two metrics, g^ v and f^, which describe two interacting 
gravitons, one of which is massive and the other one is massless. The kinetic term of 
each metric is chosen to be of the standard Einstein-Hilbert form, while the interaction 
between them is described by a local potential term U[g,f] which does not contain 
derivatives and which is expressed by a scalar function of the tensor 



7t = Vg^fa,. (2.1) 

Here g^ u is the inverse of g^ v and the square root is understood in the matrix sense, i.e. 

(7 2 r„ = 1^ = sTfau- (2.2) 

Notice that the matrix square root 7^ is well-defined when g^ v and f^ u are close to each 
other so that g' ia f au is close to the unit matrix, in which case 7^ is also close to the 
unit matrix. However, 7^ might cease to be well-defined (and/or real) for generic g,f. 

Assuming a g-matter and an f-matter interacting, respectively, only with g^ v and 
with f^ u , the action is (with the metric signature — h ++) 

%, /, matter] = — d 4 x ^g R(g) + -^ d 4 x \f^f Tl{f) 

~ j 'd 4 x ^Tg-U\g, f] + Sf [g, g-matter] + sf ] [f, f-matter] , (2.3) 

where R and 1Z are the Ricci scalars for g^ v and f^ u , respectively, k 2 = 8nG and 
kJ = 8nQ are the corresponding gravitational couplings, while k 2 = k 2 + kJ and m is 
the graviton mass. The interaction between the two metrics is given by 

4 
U = Y, h *U k (i), (2.4) 

71=0 

where b k are parameters, while £4(7) are defined by the relations 

W (7) = 1, (2-5) 

Wi(7)=£A a =[ 7 ], 



w 2 ( 7 )= E A ^ = ^(M 2 -M)' 

A<B 

W 3 (7)= E A A A B A c = i([ 7 ] 3 -3[ 7 ][7 2 ]+2[7 3 ]), 



3! 

A<B<C 

lUin) = AoA!A 2 A 3 = ^([ 7 ] 4 - 6[ 7 ] 2 [7 2 ] + 8[ 7 ][7 3 ] + 3[ 7 2 ] 2 - 6[ 7 4 ]) . 

Here A^ (A = 0, 1, 2, 3) are the eigenvalues of 7^, and, using the hat to denote matrices, 
one has defined [7] = tr(-y) = 7^ , [7*] = tr(7 fc ) = (7*)^. The (real) parameters b^ could 
be arbitrary, however, if one requires flat space to be a solution of the theory, and m to 
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be the Fierz-Pauli mass of the graviton [3J, then the five b^s are expressed in terms of 
two free parameters c 3 , c 4 as follows: 

b = 4c 3 + c 4 - 6, &i = 3 — 3c 3 - c 4 , b 2 = 2c 3 + c 4 - 1, 

&3=-(c 3 + c 4 ), 6 4 = c 4 . (2.6) 

The theory f 1 2 . 3 1) propagates 7=5+2 degrees of freedom corresponding to the 
polarizations of two gravitons, one massive and one massless. Before this theory was 
discovered [2], other bigravity models, sometimes called f-g theories, had been considered 
[30] . Such theories propagate 7+1 degrees of freedom, the additional one being the BD 
ghost [5]. 

2.1. Field equations 

Let us introduce the mixing angle rj such that n g = ncost], kj = nsmrj. Varying the 
action ( 12. 3p gives the field equations 

G» = m 2 cos 2 rj T% + T [m] » , (2.7) 

Q» = m 2 sin 2 V Tt + T [n f, (2.8) 



where G% and Q% are the Einstein tensors for g^ v and f^ u and where the two gravitational 
couplings are included into the definition of the matter sources T v and T J 1 . The 
"graviton" energy-momentum tensors obtained by varying the interaction U are 



T% = r%-U #, U = -^ rt , (2.9) 

where 

rt = {6i Uo + b 2 U 1 + b 3 U 2 + 6 4 W 3 bt 
- {6 2 Wo + & 3 W 1 + 6 4 W 2 }(7 2 )^ 
+ {6 3 W + 6 4 W 1 }( 7 3 n 

-& 4 Wo(7 4 n (2-10) 

with Wfc = Uki'j)- In deriving these expressions one uses the relations 

^5 = -f u«m a „ = -\ u«{i n )\ , (2.11) 

which can be obtained by varying the definition of 7^ and using the properties of the 
trace. The matter sources are conserved due to the diffeomorphism invariance of the 

(9) r m i (/) r m i (9) (/) 

matter terms in the action, V^ T v = 0, V/* / v = 0, where V and V are covariant 
derivatives with respect to g^ v and f^ v . The Bianchi identities for (12. 7p imply that the 

(9) 
graviton energy-momentum tensor is also conserved, V/^t = 0- Similarly, the Bianchi 

(/) 
identities imply that V/z it = 0, but in fact this condition is not independent and 
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follows from V/i T^, = in view of the diffeomorphism invariance of the interaction 
term in the action. 

If r\ — > and sin 2 r\T% — > 0, then equations (J2.8P for the f-metric decouple and 
their solution enters the g-equations ( 12. 7ft as the fixed reference metric. The g-equations 
describe in this case a massive gravity theory. For example, if f becomes flat for 77 — >• 0, 
then one recovers the dRGT theory [lj. Therefore, theories of massive gravity are 
contained in the bigravity. However, one has to emphasize that the existence of the 
massive gravity limit is not guaranteed [31] , since sin 2 rj T^ does not necessarily vanish 
when 77 — >- 0, and even if it does, the f-metric does not necessarily reduce to something 
fixed. For example, in the absence of the f-matter, f-metric becomes Ricci flat in the 
limit, but this does not guarantee that it is flat. Therefore, although the massive gravity 
can be embedded into bigravity, it is not a limit of the latter but a different theory with 
different properties. As a result, the recent observation that the massive gravity theory 
could be acausal (32] does not necessarily apply to the bigravity theory 



3. Proportional backgrounds 

The simplest solutions of the bigravity equations are obtained by assuming the two 
metrics to be proportional, 

Uu = C 2 g^ u . (3.1) 

One obtains in this case 7^ = C8^ v and so 

r"„ = (61 + 3b 2 C + 36 3 C 2 + b 4 C 3 )C5^ , (3.2) 

which gives the energy-momentum tensors 

Tt = -A g (C)S^, T» = -A f (C)5^, (3.3) 

with 

A s (C) = m 2 cos 2 rj {b + 3b ± C + 3b 2 C 2 + 63 C 3 ) , 

. 9 

sin Ti 
A f (C) = m 2 -^J- (bx + 3b 2 C + 3b 3 C 2 + b 4 C 3 ) . (3.4) 

Since the energy- momentum tensors should be conserved, it follows that C is a constant. 
As a result, one finds two sets of Einstein equations, 

g; + A g (c)s; = n , g; + A f (c)s; = i\ . (3.5) 

Since one has Q v = G^/C 2 , it follows that A/ = A g /C 2 , which gives an algebraic 
equation for C, with x — tan 2 rj, 

(b + 36i C + 3b 2 C 2 + h C 3 ) = ^ (61 + 3b 2 C + 3b 3 C 2 + WC 3 ) . (3.6) 

It follows also that the matter sources should be fine-tuned such that T ' ^ = 
T lmjp 'u /C 2 . Therefore, the independent equations are the same as in GR, 

g; + a 9 (c)s; = t^. (3.7) 
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If the parameters b^ are chosen according to Eq. (12.61) . then Eq. ()3.6p factorizes, 
= (C - 1)[(C3 + c 4 )C 3 + (3 - 5c 3 + (x - 2)c 4 )C 2 
+ ((4 - 3 X )c 3 + (1 - 2 X )c 4 - 6)C + (3c 3 + c 4 - l)x], (3.8) 



while 



V " - (1 - C)((c 3 + C 4 )C 2 + (3 - 5c 3 - 2c 4 )C + 4c 3 + c 4 - 6). 



m 2 cos 2 ?? 



Depending on values of c 3 ,c 4 , 77, the equation (J3.8P can have up to four real roots. For 
example, for c 3 = 1, c 4 = 0.3, 77 = 1 the roots are 

C = {-2.24; 1; 0.06; 3.61}, -^M = {10.12; 0; -0.50; -4.50}. (3.9) 

m 2 

As a result, there can be solutions with four different values of the cosmological constant, 
which can be positive, negative, or zero. If C — 1 then A g = A/ = and the two metrics 
and the matter sources are identical, g^ v = f^, T u = T V Setting T u = 0, the 
vacuum GR is recovered, and, in particular, flat space, g^ = f^ = 77^. 

Let us consider small fluctuations around flat space, g^ v = T]^ v + Sg^ u , f^ u = 
Vfiv + fifnv Inserting into the general equations (J2.7p . (l2.8p and linearizing with respect 
to the fluctuations, one finds that the linear combinations 

hffi = cos 77 8g^ v + sin 77 Sf^, h$ = cos 77 8f^ u - sin 77 Sg^ (3. 10) 

fulfill the Fierz-Pauli equations, 

□/»H + . . . = m \h^ - h^ V , u ), Dfcjg + . . . = 0. (3.11) 

Therefore, one can identify nffi and h^l with the massive and massless graviton fields, 
respectively. It seems that one can explicitly identify in this way the massive and 
massless degrees of freedom only for fluctuations around proportional backgrounds [34J . 
Moreover, only for proportional backgrounds the null energy condition is fulfilled both 
in the g and f sectors [35J. Apart from flat space, other solutions with proportional 
metrics do not admit the massive gravity limit with flat f-metric. 

4. FLRW cosmologies with non-bidiagonal metrics 

Let us now make a symmetry assumption and choose both metrics to be invariant under 
spatial SO (3) rotations. Since the theory is invariant under diffeomorphisms, one can 
choose the spacetime coordinates such that the g-metric is diagonal. However, the 
f-metric will in general contain an off-diagonal term, so that 

ds 2 g = _ Q2 rft 2 + ^2^2 + ^2^2 ? 

ds) = - (aQdt + cNdrf + (cQdt - bNdr) 2 + u 2 R 2 dtt 2 . (4.1) 

Here Q,N,R, a, b, c, u depend on t, r and dVL 2 = d'd 2 + sin 2 'ddtp 2 . One can take 

( a cN/Q \ 

-cQ/N b 

u 

\ u ) 



7t = \/9 m fa V 



(4.2) 
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whose eigenvalues are 

A ,i= Ua + b± v /{a-by-AcA, X 2 = A 3 = u. (4.3) 

Inserting this to (J2.5P gives 

U 1 = a + b + 2u, U 2 = u{u + 2a + 26) + ab + c 2 , 

U 3 = u(au + bu + 2ab + 2c 2 ), U A = u 2 (ab + c 2 ) . (4.4) 

Notice that, although the eigenvalues (14.31) can be complex- valued, the £4's are always 

real. It is now straightforward to compute the energy-momentum tensors T^ u and T^ 

defined by (I2.9IUI2.10I) . In particular, one finds 

cN 
T Q r = C —[b 1 + 2b 2 u + b,u 2 }. (4.5) 

It will be assumed in what follows that the g-metric is either static or of the FLRW type, 
in which cases there is no radial energy flux and T° r should be zero. Therefore, either 
c should vanish, or the expression in brackets in (14. 5B vanishes. The former option will 
be considered in the next section, while presently let us assume that c^ and 

&i + 2b 2 u + b 3 u 2 = 0. (4.6) 

This yields 



u = i (-b 2 ± y/b* - hhj . (4.7) 

Notice that u was a priori a function of t,r, but now it is restricted to be a constant. 
Using this, one finds that T° = T r r = —X g and T° = T\ — —A/ where 

\ a j_ oh j_ u 2 x & 2 + 26 3 « + b 4 u 2 

\q = D + 2biu + b 2 u, Xf = 5 . (4.8) 

u l 

(a) 
The conditions Vp 2^ = reduce in this case to the requirement that T° — T\ should 

vanish. On the other hand, one finds 

T° - T% = (b 2 + b 3 u)[(u - a)(u -b) + c 2 ], (4.9) 

and since this has to vanish, either the first or the second factor on the right should be 
zero. The former case was considered in [H],[23] (see also [19]). However, requiring that 
b 2 + b 3 u = constrains the possible values of the theory parameters b& hence solutions 
obtained in this way are not general. Let us therefore require that J2TJ [22] 

(u-a)(u-b)+c 2 = 0. (4.10) 

Note that in this constraint u is a constant, but a, b and c are still (unknown) functions 
of t,r. In view of this, one has T° = T^ and 1~° = T\ hence both energy-momentum 
tensors are proportional to the unit tensor, T^ = — X g S^ and T^ = — Xf8%. The field 
equations ( 12. 7p then reduce to 

G{ + A g 6? x = T^ )p x , 

Q{ + A A P = r M A p , (4.11) 
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with K g = m 2 cos 2 7/ A and Aj = m 2 sin 2 r\ Xf. As a result, the two metrics seemingly 
decouple one from the other, and the graviton mass gives rise to two cosmological terms. 
Unlike in the case of proportional metrics, no fine tuning between the two matter sources 
is needed. However, one has to remember that solutions of (14. lip should fulfill the 
consistency condition ( I4.10p . If the parameters bk are chosen according to (12. 6p then 
X g + u 2 X f = —{u - l) 2 , therefore, if A g > then A f < 0. 

Let us consider time-dependent solutions of the FLRW type and assume the g- 
matter to be a perfect fluid, T ^ = diag[— p g (t),P g (t),P g {t), P g {t)\, while the f-sector 
can be chosen to be empty, T A p = 0. Then the solution of (14.1 ip is chosen to be 

ds 2 g = -dt 2 + a 2 (£) f ^ 2 + r 2 dtt 2 \ , k = 0, ±1, 

ds 2 = - A(U) dT 2 + 4^r + U 2 dn 2 , A = 1 - ^ U 2 , (4.12) 

J A(U) 3 

where a fulfills the Friedmann equation, 
„ a 2 + k 



a 2 



A g + p. (4.13) 



The g-metric describes an expanding FLRW universe containing the matter p and the 
positive cosmological term A g . It can be of any spatial type - open, closed or flat. At 
early times, when a is small, the matter term p is dominant, while later A g dominates 
and the universe enters the acceleration phase. The f-metric is the anti-de Sitter one 
expressed in static coordinates. One has Aj ~ sin 2 r\ — >■ when r\ — V 0, hence the 
f-metric becomes flat in this limit. Therefore, the solutions apply both in the bigravity 
theory and in dRGT massive gravity. 

4-1. Imposing the consistency condition 

The described above effective decoupling of the two metrics has been observed by several 
authors, in the massive gravity theory [151 021 HH1 [20] and also in the bigravity [23J. 
However, such a decoupling is only possible if the consistency condition (I4.10}1 is fulfilled, 
which requires solving a complicated PDE (see Eq. (l4.14p below), and this has not always 
been done. In [15] the solution was obtained for p g = when the g-metric is the de 
Sitter one expressed in static coordinates, in which case the PDE becomes an ODE. A 
solution for a more general FLRW g-metric for k = and for special values of bk was 
found in [18], while the general case was considered in [2"T] |2"2"]. 

Let us notice that the two metrics in (14.121) are expressed in two different coordinate 
systems, t, r and T, U, whose relation to each other is not yet known. One has T = T(t, r) 
and U = U(t, r), so that dT = Tdt + T'dr and dU = Udt + U' 'dr. Inserting this into 
the f-metric in (14.121) and comparing with the f-metric from (14.11) . one finds the metric 
coefficients a, b, c in ( 14. ip expressed in terms of the partial derivatives of T, U. One also 
finds that U = uR = uar. Inserting the result into (I4.10p gives a non-linear PDE, 

A [ a Vl - kr 2 (UT' - TU') - w 2 a 2 l = u 2 a 2 A + A_ , (4. 14) 
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with A± = a (AT ±U) + y/l - kr 2 (W ± AT'). Since u, a, U, A(U) are already known, 
this equation determines T(t,r). 

Let us consider the rj — y limit, when A/ = and A = 1. Then one can find exact 
solutions of (I4.14p . If fc = then 

T{t,r) = qj j + (^ + qAa, (4.15) 



where q is an integration constant. This solution agrees with the one obtained in [18 
for C3 = C4 = 0, u = 3/2. For k = ±1 one has 



T(£, r) = vV + ku 2 / Va 2 + k dt + gaVl - kr 2 . 



(4.16) 



If A/ 7^ and A^l then exact solutions of (14. 14ft are unknown, however, at least when 
Af is small, solutions can be constructed perturbatively as T = T + ^ n>1 (— Af/3) n T n . 
Here To corresponds to zero order expressions ( 14.151) , ( 14.161) , while the corrections T n can 
be obtained by separating the variables with the ansatz T n = Y^=o fm(t)(l — kr 2 ) m ^ 2 . 

This completes the construction, since all field equations and the consistency 
condition are fulfilled. Summarizing, one obtains self-accelerating FLRW solutions of 
all spatial types, equally valid in the bigravity and massive gravity theories. In the 
latter case this exhausts all possible homogeneous and isotropic cosmologies. Somewhat 
confusingly, these solutions are sometimes called in the literature inhomogeneous, 
because the two metrics are non-bidiagonal and T, U, which play the role of Stiickelberg 
scalars in the massive gravity limit, are inhomogeneous, as they both depend on r. 
It is then expected that fluctuations around the FLRW backgrounds should show an 
inhomogeneous spectrum, although this effect will be proportional to m 2 and so will be 
small [IE]. 

In the spatially open case, k = —1, choosing in (14. 16ft q — u, yields T = uay/l + r 2 . 
Inserting this to ( I4.12p . together with U = uar, gives the flat f-metric which turns out 
to be diagonal both in the T, U and t, r coordinates, 

ds 2 = -dT 2 + dU 2 + U 2 dn 2 = u 2 a 2 (-% dt 2 + -^- + r 2 dn 2 ) , (4.17) 
j \ a 2 I + r 2 J 

so that the g and f metrics happen to be bidiagonal in this case. This particular solution 
was found in [3S] • It is unclear if it should be physically distinguished when compared 
to the other solutions in (I4.15l) . ()4.16p for generic k,q. The analysis of [37] shows that 
this solution is unstable at the nonlinear level, while the case of generic k, q has not yet 
been studied. 

5. FLRW cosmologies with bidiagonal metrics 

All solutions considered above are described by the same metrics as in GR, and 
the graviton mass manifests itself only as the effective cosmological term(s). More 
general solutions are obtained assuming that both metrics are simultaneously diagonal 
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dr 2 



rf S 2 = _ dt 2 + a(t) 2 ^^_ + r 2 dfi 2 



dr 

The Gg and ^ equations (EZD,(E8D then read 

a 2 _ A g (Q+p g _fc_ ^b 2 A f (0+pf k_ 
a 2 3 a 2 ' „4 2 b 2 3 b 2 

where £ = b/a and where A 9 , A^ are defined in (j3.4p . The V^t = conditions reduce 
to 



rfs 2 = - ^(t) 2 rft 2 + b(t) 2 + rW . (5.1) 



(5.2) 



b-^la 



5.1. Generic solutions 



(p 1 + 2 b 2 i + b z e) =0. (5.3) 



Let us assume that the first factor in (15.31) vanishes, b — Ak = 0. Comparing this to 
(15. 2p gives the algebraic relation A g (£) + p g = £ 2 (Ay(£) + pf), or, setting for simplicity 

|e 4 +(&2-f&4)£ 3 + (&l-^3)£ 2 + (|-X&2 + ^V)£=^,(5.4) 
o o o o COS Tf o 

with x — tan 2 rp Assuming the equation of state P g = w g p g one has p g = p g sr 3( - 1+W3 \ 
and so that Eq. fl5.4j) gives an algebraic relation between a and £, whose solution is £(a). 
Inserting £(a) into the first equation in (15.21) gives the Friedmann equation, 

a 2 + V(a) = -k , V(a) = -y (A 3 (0 + p g ), (5.5) 

where £ = £(a) and p g = p 9 (a), hence a is confined to the region where V^(a) < —k. 

Depending on choice of the solution of (j5.4p . the shape of the potential V(a) can be 
different. The simplest solutions are obtained for b\ = 63 = [21]. Making the parameter 
choice (12.61) . the self-accelerating solutions are found only if b\ — 3 — 3c 3 — c 4 < [23J. 
In this case Eq. fj5.4p has only two real roots £i(a) and ^2( a ) which exist for any positive 
a, A g (£i) approaching a positive value at large a while A g (£ 2 ) being negative. This 
corresponds, respectively, to the accelerating and recollapsing solutions (see FigJT]). For 
61 > there could be more complex solutions for which the universe size oscillates around 
a finite value [23]. When Eg. (15.51) is solved, the f-metric coefficients are b = a£(a) and 
A = b/a. 

The perturbative analysis reveals that the self-accelerating solutions are stable 
[39| HOI HI] , hence they can be used for interpreting observational data 



5.2. Special solutions 

Let us now assume that the second factor in (15.31) vanishes, 61 + 26 2 £ + b^ 2 = 0. Notice 
that this is the same equation as in (14.61) . It follows that £ = b/a is constant, therefore 
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Figure 1. The two solutions £1, £2 of (J5.4J) and the corresponding effective 
potentials V(a) for the same parameter values as in (|3.9p . One has V% < — 1, 
hence for any fe = 0,±1 solutions of (|5.5j) travel over the Vi potential barrier 
and enter the accelerated phase for large a. Since A 9 (^2) < 0, Vi — > +00 for 
large a and solutions of (|5.5p bounce from the potential barrier and recollapse. 



a/a = b/b. Equations then again reduce to (15. 5p . where A g is now constant. Combining 
the two equations ( 15. 2p gives 

£ 2 a 2 



A 1 



(5.6) 



£ 2 a 2 (A / + P/ )/3-fc' 

and this expression can be positive only in special cases. In the massive gravity limit, 
where pf = r\ = Af = 0, A 2 is positive only for k = — 1, which reproduces the solution 
(I4.17p . and A does not exist for k = 0,1, hence there are no FLRW cosmologies with 
bidiagonal metrics in these cases |18| . 



6. Anisotropic cosmologies 



Solutions considered in the previous section can be generalized to the anisotropic case 
by choosing the two metrics to be bidiagonal and of the same Bianchi type [27] , 

ds 2 g = -a 2 dt 2 + h ah u a ®u b , ds 2 f = -A 2 dt 2 + n ab u a (g)u b . (6.1) 

Here a, h a b, A, % a b depend on t, and the vectors e a dual to the one-forms u a generate a 
three-parameter translation group acting on the 3-space, [e a , e&] = C c ab e c . The structure 
coefficients are C c ab = n cd ed a b with n ah = diag[n^, n^ 2 ',n^ 3 '], where nS a ' assume the 
values shown in Table 1. 

This choice of the structure coefficients corresponds to the Bianchi types of class A, 
in which case one can choose h a b = diag[af , af, af] and 7i a b = dia.g[A 2 ,A 2 ,A 2 }, and 
this guarantees that G° = Q® = 0, so that radial energy fluxes are absent. Choosing 



[ai,a 2 ,a 3 J 
[Ai,A 2 ,A3\ 



e nr e /3++%/3/3_ e P+-V3P- e -2/3+l 



Wr e B++V3B- e B+-V3B- 



-2B 4 



the spatial 3-curvature of the g-metric is 



(3)_ 2n^nW 

R — n 



2aja2 a : 



(n^al + n^al-n^at)' 



(6.2) 



(6.3) 
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I 


II 


VIo 


VIIo 


VIII 


IX 


n w 





1 


1 


1 


1 


1 


n (2) 








-1 


1 


1 


1 


n (3) 














-1 


1 



Table 1. Values of n^ for the Bianchi class A types 



(3) 
while the 3-curvature 1Z of the f-metric is obtained from this by replacing a a —¥ A a and 

The potential (j2.4p is U^/^g = aU g + AUj where 

Ug = boe ™ + he w + he w +2 n ^2(B + - M + 2e B + -(s + cosh[v /3 ( #_ _ /3 _ ) ]^ 

+ b 2 e 2W+n [e 2 ^-^ + 2e-^-^ cosh[V3(B_ - /?_)]) , (6.4) 

while Uf is obtained from this by replacing && — > bk+x- The equations for the g-metric 
read 




fa 



+ 



a 



P- 



a 



,-3fii 



(3) 



2cos^e Ml U g - R +2p g 

(3) 



(r\T t 
— + 3aU g ) -2ae su R +3ae 3U {p g - P g) 



l__d_ 



9 ,0 (3) 

2 cos 2 rjU -a e Mi R 



(6.5) 



while the equations for the f-metric are obtained from these by replacing Q — > W, 

(3) (3) (g) 

f3± — > W±, R—^TZ, U g —¥ Uf and cost] -> sin 77. The conditions V/^t = reduce to 



n ( i7^+B +7 ^- + s_^£-]r„ 



aw 



9/3, 



«9/3_ 



^'^^^^ 



W 



/■ 



All quantities in these equations are dimensionless, assuming 1/m to be the length scale, 
while the energy density is measured in units of m 2 M 2 [ . 

For the Bianchi types I or IX one can set anisotropies to zero, (3± = B± = 0, and 
then, with a = 2e n and b = 2e w , the above equations reduce to the FLRW equations 
(I5.2p .f l5\3|) for k — 0, 1 (the spatially open k — — 1 FLRW case is contained in the Bianchi 
V type, which does not belong to the class A). 

In the Bianchi I case it is consistent to set /3± = B± and then the two metrics 
are proportional as described in section [3j f^ v = C 2 g^ u , if only the matter sources are 
fine-tuned such that pf = p g /C 2 . One has 



dst 



. di * + e ™ ^ + +V3P- dx 2 + e 2? + -V3P- dx 2 + e -m dx fj 



(6.6) 
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where Q, (3± fulfill (assuming that a = 1) 



Q 2 



[<r. 



a 2 )e- 6n 



+ - (A 9 (C) + p 9 ) 



/3± = cT ± e~ Ml , (6.7) 

with constant a±. Here C, A P (C) = 3H 2 are defined by (El} , ([3J| . At late 
times solutions of Eqs. f)6.7p approach configurations with constant anisotropies, Q = 
Ht + 0(e~ 3Ht ), P ± = /3±(oo) + 0(e~ 3Ht ). 

It turns out that solutions for all other Bianchi types and for non fine-tuned sources 
also evolve into configurations with equal and nonvanishing anisotropies, (3± = B±. This 
can be seen by applying a numerical procedure [27], whose input configuration at the 
initial time moment t = is an anisotropic deformation of a finite size FLRW universe. 
The initial value of Q is zero, hence the initial universe size e n ~ 1 (in 1/m units), while 
the initial anisotropies are (3±, B±, (3±, B± ~ 1CT 2 . The initial values of Q, W, W are 
determined by resolving the first order constraints contained in the field equations. It 
is assumed that the f-sector is empty, pf = 0, whereas the g-sector contains a radiation 
and a non-relativistic matter, p g = 0.25 x e~ AQ + 0.25 x e~ 3f \ so that the dimensionful 
energy w?M 2 x p g ~ 10~ 10 (eV) 4 , assuming that m ~ 10 _33 eV. 

The numerical extension to t > of the t = initial data reveals that, for all 
Bianchi types, the universe approaches the state with proportional metrics in which 
the expansion rate Q is constant, H, and the anisotropies f3± = B± are also constant. 
For the Bianchi I solutions the constant anisotropies can be absorbed by redefining 
the spatial coordinates, but not for other Bianchi types. This means that the universe 
generically runs into anisotropic states at late times. These properties of the solutions 





0.5 



Figure 2. The Hubble parameter Q for all Bianchi types of class A (left) and 
the anisotropy amplitude for the Bianchi IX solution (right). 



are illustrated in Figf2J The relative shear contribution to the total energy density 

S = 



$+ + /3-/^ tends to zero at late times, however, if only just one or two Hubble 
times have elapsed since the beginning of the current phase of the universe acceleration, 
E should not necessarily be small. It turns out that at late times the anisotropies 
oscillate around their asymptotic values. Linearizing the field equations with respect to 
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small deviations form the proportional background f^ u = C 2 g^ u to which the solutions 
approach, one finds that 

$± ~ B± ~ e -3Ht/2 cos(Hut), (6.8) 

where u; can be expressed in terms of C,bj,,r) [27] . It follows that the shear energy 
is /3i + P 2 _ ~ e~ 3n ~ 1/a 3 , while in GR it decreases as ~ 1/a 6 . The fact that the 
shear energy in the bigravity theory shows the same fall-off rate as a cold dark matter 
suggests that the latter could in fact be an anisotropy effect, although it is unclear if 
this interpretation can also explain the dark matter clustering. 







ln(cti^-^__ 


---^ ^^^"T/ 7 " 
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^^J^^^C^/ - 
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/ 
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-40 -30 -20 -10 

Q. 



Figure 3. Near singularity behavior for the Bianchi I (left) and Bianchi IX 
(right) solutions. 

When continued to the past, the solutions hit a singularity where both e Q and e w 
vanish (see Figj3]). The Bianchi IX solutions, when approaching singularity, show the 
typical billiard behavior characterized by a sequence of Kasner-type periods [13]. During 



each period one has a p oc t Pa where p\ + p 2 + P3 



Pi + v\ + vl 



1, such that two of 



the three a a 's grow in time whereas the third one decreases, until the next period starts 
when one of the growing amplitudes becomes decreasing (see Figj3]). 

7. Black holes, lumps, stars and the Vainshtein mechanism 



Historically, the first bigravity black holes were obtained in the generic f-g bigravity 
theory [30] for non-bidiagonal metrics |44| . Their counterparts in the ghost-free bigravity 
under consideration are described by the ansatz (14.11) . The field equations then reduce 
to ( 14. lip (without matter terms), and the solution is the Schwarzschild-de Sitter metric, 



dst 



M 2 + ^ +r 2 do2 



2M 



A g r 2 



D ■■-' ° = 1 ---^ «.±) 

with A g defined after Eg. (14. 111) . The f- metric is still expressed in terms of T, U by (I4.12p . 
and to fulfill the consistency condition (14.101) the procedure is the same as in Section HI 
which gives U = ur and T = ut — u J D ^^ dr with u from ( 14. 7p [28J. Since the f- metric 
becomes flat for r\ — y 0, this solution describes black holes also in the dRGT massive 
gravity [13 US] (see also [UJ ) . 
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These solutions and their generalizations for a nonzero electric charge [H] exhaust 
all known black holes in the massive gravity theory. In particular, there are no 
asymptotically flat black holes in this case (see [46] for a recent discussion). However, 
in the bigravity one finds more solutions when the metrics are bidiagonal [28] . 



ds 2 = Q 2 dt 2 - d C- rW, ds 2 = A 2 dt 2 - ^-dr 2 - U 2 dtf . (7.2) 

y iv Y 

Here Q, N, Y, U, a are 5 functions of r which fulfill the equations 
G° = m 2 cos 2 7] T °, G r r = m 2 cos 2 rj T r r , 

Qq = m 2 sin 2 rjT^, Q r r = m 2 sin 2 77 77, 



->r/ 1 ^ trrrr rpO\ 1 ^ (rrfd 



T; + ^ (T; " T ° } + r W " T;) = °- (7 - 3) 

The simplest solutions are again obtained for the proportional metrics, f^ u = C 2 g^ v , 
with g^ u given by ( 17. ip where A g = A g (C) is defined by (I3.4p .( )3l6|) . Since A g can 
be positive, negative, or zero, there are the Schwarzschild (S), Schwarzschild-de Sitter 
(SdS), and Schwarzschild- ant i-de Sitter (SAdS) black holes. Let us call them background 
black holes. More general solutions are obtained by numerically integrating Eqs. (l7.3p . 

In doing this, it is assumed that the g- metric has a regular event horizon at r = r^, 
where Tjf and the curvature are finite. A detailed analysis shows [28] that the horizon 
should be common for both metrics, therefore amplitudes Q 2 , N 2 , A 2 , Y 2 should have a 
simple zero, while U, U' should be non-zero at r = r^. Next, one finds that the boundary 
conditions at the horizon comprise a one-parameter family labeled by u = [/(r^/rh, 
the ratio of the event horizon radius measured by the f-metric to that measured by 
the g-metric. Finally, it turns out that the horizon surface gravities and temperatures 
determined with respect to both metrics are the same [18] . 

Choosing a value of u and integrating numerically the equations from r = Th towards 
large r, the result is as follows [28J. If u = C where C fulfills Eq. (l3.6p then the solution 
is one of the background black holes described above. If u = C + 5u with a small 5u, 
then one can expect the solution to be the background black hole slightly deformed by a 
massive graviton 'hair' localized in the horizon vicinity. This is confirmed for the SAdS 
solutions (A g < 0), which indeed support a short massive hair and show deviations from 
the pure SAdS in the horizon vicinity, but far away from the horizon the deviations tend 
to zero (see Fig|4]). 

However, the argument does not work for A g > 0. When one deforms the S 
background by setting u = r^ + 8u for however small 5u, the solutions first stay 
very close to Schwarzschild (until ln(r/r/j) ~ 4 in FigH|). However, for large r they 
deviate from the background and show a completely different asymptotic behavior at 
infinity, characterized by a quasi- AdS g-metric and a compact f-metric giving a finite 
value for the volume of the 3-space [28] . Therefore, the only asymptotically flat black 
hole in the bigravity theory is pure Schwarzschild, and deforming it results in loosing 
the asymptotic flatness. Similarly, trying to deform the SdS background produces a 
curvature singularity at a finite proper distance away from the black hole horizon, hence 
the only asymptotically de Sitter black hole is the pure SdS. 
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Figure 4. Hairy deformations of the SAdS background, where Aq,Nq,Qq,Yq 
correspond to the undeformed solution (left); and of the S background (right). 



In the shrinking horizon limit, r^, — > 0, the black hole 'hair' does not disappear 
but becomes a static 'lump' made of massive field modes. Such lumps are described by 
globally regular solutions for which the event horizon is replaced by the regular center 
at r = 0, while at infinity the asymptotic behavior is the same as for the black holes 
[28] . None of the lumps are asymptotically flat (apart from the flat space). 

Asymptotically flat solutions can be obtained by adding matter. Suppose that 
the f-sector is empty, while the g-sector contains T [m '^ u = diag[— p(r), P(r), P(r), P(r)] 
with p = p^Oi^r* — r), corresponding to a 'star' with a constant density p* and a radius 
r*. Adding this source to the field equations (I7.3I) . assuming the regular center at 
r = and integrating towards large r, one finds a solution for which both metrics 
approach Minkowski metric at infinity. Introducing the mass functions M g ,Mf via 
g " = N 2 = 1 - 2M g (r)/r and f rr = Y 2 /U' 2 = 1 - 2M / (r)/r, one finds that M gi M f 
rapidly increase inside the star, while outside they approach the same asymptotic value 
M g (oo) = Mj (oo) ~ sin 2 r] (see Figj5]). For r] = n/2 the g- metric is coupled only to 
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Figure 5. Profiles of the asymptotically flat star solution sourced by a regular 
matter distribution. 



the matter and is described by the Schwarzschild solution, M g (r) = p*r 3 /6 for r < r* 
and M g (r) = p*rf/6 = Madm for r > r*. For r\ < n/2 the star mass Madm is partially 
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screened by the negative graviton energy. For rj = the f-metric becomes flat, so that 
Mf = 0, while M g asymptotically approaches zero and the star mass is totally screened, 
because the massless graviton decouples and there could be no 1/r terms in the metric. 

If the graviton mass is very small, then the m 2 T^ v contribution to the equations 
is small as compared to T „, and for this reason M g rests approximately constant for 
r+ < r < ry ~ (Madm/^ 2 ) 1 ^ 3 - This illustrates the Vainshtein mechanism of recovery of 
General Relativity in a finite region [29] . This mechanism has also been independently 
confirmed within the generic massive gravity theory with the BD ghost [4*9] 150] , and in 
the dRGT theory [51] . It is worth noting that the mechanism works only in the presence 
of a regular matter source, and not for vacuum black holes or lumps. 

To recapitulate, globally regular and asymptotically flat stars exist both in the 
bigravity theory and in the dRGT massive gravity. In the latter case there are also SdS 
black holes, but there are no asymptotically flat black holes. There are hairy black holes 
in the bigravity theory, but the only one which is asymptotically flat is the standard 
"hairless" Schwarzschild black hole. 

8. Concluding remarks 

Summarizing what was said above, all known cosmologies and black holes in the 
dRGT massive gravity theory are obtained within the inhomogeneous Stiickelberg field 
approach, assuming non-bidiagonal metrics [36J. Perturbations of the cosmological 
solutions were studied in [521 |53j [511 EH] ■ O ne could perhaps also mention a peculiar 
FLRW solution [17] with a degenerate f-metric, although its g-metric is regular and is 
sourced by the effective T° = ^ n=0 a n a_n with constant a n (see also [56| 157]). 

Although massive gravity solutions with bidiagonal metrics similar to (14.171) are 
very special and hence are unlikely to be physically important, they are popular due to 
their simplicity. They have been studied also in the modified versions of the theory, as 
for example in the extended massive gravity [581 159] 160] where the graviton mass m is 
promoted to a dynamical field [18] [61] , or in the massive gravity theory with a (quasi) 
dilaton field J62J [63]. Choosing the fixed f-metric to be non-flat, as for example de Sitter 
or FLRW, also leads to new cosmological [Ml ESI [661 E3 EH [69] and black holes solutions 
[70] . and allows one to study [7_H [721 E3 the Higuchi bound [7_I], that is the limit where 
the scalar polarization of the massive graviton decouples. In addition, 0(4)-symmetric 
instanton solutions in the Euclidean version of the massive gravity theories have been 
considered as well [7_5l [76] . 

Changing the reference metric from case to case does not seem to be natural. If 
one wants to consider different possibilities for f, then it is logical to make it dynamical, 
which leads to the bigravity. The massive gravities with a fixed f then can be recovered 
by choosing a source for f, for example a constant pf = Pf > 0, and taking the limit 
T] — > 0. For rj y^ both metrics are dynamical, which gives rise to a much richer solution 
structure than in the massive gravity theories. 
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